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This paper aims to use a new technique of computed tomography (CT) scan 
image processing to correlate the image analysis with sinonasal symptoms. 
A retrospective cross-sectional study is conducted by analyzing the digital 
records of 50 patients who attended the ear, nose and throat (ENT) clinics at 
King Abdullah University Hospital, Jordan. The coronal plane CT scans are 
analyzed using our developed software. The purposes of this software are to 
calculate the surface area of the nasal passage at three different levels visible 
on coronal plane CT scans: i) the head of the inferior turbinate, ii) the head 
of the middle turbinate, and iii) the tail of the inferior turbinate. We employ 
image processing techniques to correlate the narrowing of nasal surface area 
with sinonasal symptoms. As a consequence, obstruction in the first level is 
correlated significantly with the symptoms of nasal obstruction while the 
narrowing in the second level is related to frontal headache. No other 
significant correlations are found with nasal symptoms at the third level. In 
our study, we find that image processing techniques can be very useful to 
predict the severity of common nasal symptoms and they can be used to 
suggest treatment and to follow up on the case progression. 
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1. INTRODUCTION 


The nasal cavity contains important structures, such as vestibule, nasal valve, septum, olfactory 
region, paranasal sinuses, and turbinates as shown in Figure 1, which can have different sizes and shapes in 
different people. However, in general, the main structural characteristics of the nose are consistent between 
all individuals. Simple changes in the nasal cavity’s normal anatomy may alter the airflow and lead to 
problems such as snoring, nasal obstruction and even facial pain [1]. 

The nose and paranasal sinuses are integral components of the body, which are affected by a number 
of inflammatory and neoplastic conditions affecting about 50 million individuals annually [2], [3]. 
Diagnosing issues like nasal obstruction, headache, anosmia, snoring, epistaxis, oral symptoms, facial 
swelling, orbital symptoms, and ear symptoms, is based on history and abnormal physical exams, such as 
purulent discharge, pale edematous turbinate, and nasal masses, [4], [5]. Nasal obstruction is considered a 
symptom of high subjectivity. Moreover, it is among the most widely faced complaints in otolaryngology 
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clinics. Many patients had nasal surgery to improve nasal patency and airflow. Regarding these points, 
doctors have long been interested in techniques providing objective evaluations for this complaint [6]. 


Olfactory region 


Auditory tube 
opening 


Vestibule 
Nasopharynx 


Respiratory region 


Figure 1. Nasal cavity anatomy [1] 


The radiological evaluation of sinonasal diseases is very critical as the clinical findings might be 
unclear in many cases [7], [8]. The imaging techniques available are plain radiography, computed 
tomography (CT), positron emission tomography (PET), and magnetic resonance imaging (MRI). There are 
pros/cons for each modality. For example, since plain films can only give general or high-level ideas of 
pathology and anatomy, they are not considered part of the primary imaging [9]. On the other hand, CT and 
MRI can show fine anatomic and pathological details in serial sections [9]. 

An excellent imaging technique for assessing the sinonasal cavities is CT scanning due to its 
accurate assessment of the nasal septum, the paranasal sinuses, in addition to the turbinate and olfactory 
region. For the diagnosis of many abnormalities in the sinus (for both adults and children), CT scanning is 
viewed by many as the radiographic gold standard [10]. Moreover, the CT scan is the best option in defining 
the complex anatomy and variations that are inaccessible by examination or endoscopy because of its 3D 
high resolution [11]. 

In many cases, the patients complain from various sinonasal symptoms, and, once CT scans are 
done for them, nothing abnormal is discovered regarding the sinuses. Thus, based on this point, we formulate 
our research objective to correlate between symptoms, such as rhinorrhea, hyposmia, facial pain, nasal 
obstruction and snoring, with the diameter of the nasal cavity measured in a CT scan at three different levels 
with respect to the nasal septum: (L1) the head of the inferior turbinate, (L2) the head of middle turbinate and 
finally (L3) the tail of inferior turbinate. 

The rest of this paper is organized as follows: section 2. Describes the methodology we use in 
thiswork. The results and discussion of our analysis is discussed in sections 3 and 4. Finally, a conclusion of 
the paper is presented in section 5. 


2. RESEARCH METHOD 

This cross-sectional study is conducted by analyzing the database of 50 patients who attended the 
ear, nose and throat (ENT) clinics of King Abdullah University Hospital (KAUH), Irbid, Jordan, from 
January 2018 to May 2018. The size of this dataset may seem small, but that is a common characteristic of 
medical image datasets [12], [13]. After interviewing them and asking them to fill out a questionnaire, CT 
scans were conducted for patients on the first visit to our clinic before being prescribed medical treatments. 
All patients with evidence of opacification on CT sinuses or had the previous sinus or septal operations were 
excluded from the study. To create our dataset, we have divided our work into three phases: dataset 
collection, images preprocessing, and region of interest (ROI) selection. These phases are detailed in the 
following subsections. 


2.1. Dataset collection 

This study is concerned with finding a relation between symptoms that the patient had (including 
nasal obstruction, headache, rhinorrhea, hyposmia, and snoring) and the area of nasal passage at three 
different levels as demonstrated in Figure 2 [14]. The head of the inferior turbinate represents the first level. 
The head of middle turbinate is the second one and the tail of inferior turbinate is the final one. To construct 
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this dataset, several procedures have been conducted with targeted patients. The contributing doctors 
diagnosed the patient to decide whether the case is appropriate for the study or not. Then the patient is 
questioned by the doctors to collect the required information about the patient in addition to the patients’ CT 
scans. 

The information collected from patients distributed into three groups of features. The first and 
second groups are gathered from patient interviews and asking the patients to fill out a questionnaire. The 
third group includes the information gathered from the patient CT scan (the air percentage in the paranasal 


cavity). 


F Frontal Sinus 

/ Nasal Cavity 

/ / Superior Turbinate 
A / — Middle Turbinate 
/ ;— Inferior Turbinate 


Spenoid Sinus 


Nasopharynx 
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Eustachian Tube Orifice 


Fossa of Rosenmuller 


Figure 2. Turbinate’s anatomy [14] 


The second group of features concerns about patient illness severity. This group has nine features; R 
obs (obstruction in the right side of the nose), L obs (obstruction in left side of the nose), R rhino (rhinorrhea 
in the right side of the nose), L rhino (rhinorrhea in the left side of the nose), R headache (headache in the 
right side of the head), L headache (headache in the left side of the head), Hyposmia (decreases ability to 
smell), Sneezing and finally Snoring. The values of these nine features range between 0 and 10, where the 
patient weights the severity of the symptoms that he/she has according to the visual analog scale (VAS). VAS 
is a scale or a tool that helps in assessing the symptoms intensity of the patient. This scale is a line that is 
around 10cm with anchors labeled with statements starting from the left with no pain to the right with worst 
possible, unbearable pain as shown in Figure 3. The patient is asked to decide the level of his/her pain by 
placing a mark on one of those anchor statements. 


Visual analog scale 


No Worst pain 
pain imaginable 


Verbal descriptor scale 


No Mild Moderate Severe Very Worst 
pain pain pain pain severe pain 
pain imaginable 


Numeric rating scale 


No Moderate Worst 
pain pain pain 
imaginable 


Figure 3. A visual analog scale [15] 
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The last group of attributes depends on computing the air percentage in the nasal cavity area. The 
turbinates are structures that exist in the lateral nasal wall. Clarification images of this area are shown in 
Figures 4 and 5. The air percentage in this area can be calculated from the coronal view of the patient’s CT 
scan of the area. The CT scan is usually stored as digital imaging and communications in medicine (DICOM) 
files. These DICOM files usually have tens of coronal slices of the area. We collect sinus CT scan images 
and targeted three zones for detailed elaboration. The first zone (L1) is located anteriorly between the head of 
inferior turbinate and septum medially. Zone two (L2) is in the middle between the head of the middle 
turbinate laterally and nasal septum medially. The third zone (L2) is posteriorly at the level of the posterior 
choanae with the tail of the inferior turbinate. We study images of three levels and select the narrowest zone 
for each level based on dealing with the nasal cavity as a tube and the narrowest site will be the most 
symptomatic zone. The sinus results are in the three levels as previously mentioned. This group has 6 
numeric attributes; R1 represents the percentage of air in L1 at the right side of the nasal septum, L1 
represents the air percentage in L1 at the left side of the nasal septum, R2 represents the percentage of air in 
L2 at the right side of the nasal septum, L2 represents the air percentage in L2 at the left side of the nasal 
septum, R3 represents the percentage of air in L3 at the right side of the nasal septum, and L3 represents the 
air percentage in L3 at the left side of the nasal septum. Through CT scan the three levels of turbinate can be 
observed within different slices; Figure 4 shows the three levels in the coronal plane and Figure 5 shows the 
three levels in the axial plane. 


level 3 


Figure 3. A coronal CT scan illustrating the three different levels that we consider for our correlation analysis 
of the patients' nasal symptoms 


Figure 4. An Axial computed tomographic (CT) scan illustrating the three different levels that we depend on 
the CT to correlate it with patients’ nasal symptoms 


2.2. Images preprocessing 

All of the images are obtained using CT scans and are stored as DICOM files. A DICOM file 
consists of 2D arrays with pixel intensities [15], [16]. Our developed tool takes a DICOM file as input and 
shows it as a grey-scale image. This type of files can not be viewed directly as another type of images in the 
MS Windows operating system since it consists of a lot of metadata [17]. Figure 6 shows how the user can 
load and open DICOM files using our developed tool. The figure shows the open dialog that prompts the user 
to determine the DICOM file location. This step done by clicking on the “Load image” button. The DICOM 
image loaded and displayed in the leftmost box as the figure shows. The reading process is done using the 
pydicom python package for DICOM files version 1.0.2. 

The next step is to enhance the quality of the images by removing noise using the median filter with 
kernel size equal to 3; median filter is used to remove ‘salt and pepper' noise type while saving the edges. 
This filter passes through the image pixel by pixel and replaces the gray-level value of each pixel by the 
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median of its neighbors [18], [19]. The next step is to resize images to have the same size for all of the patient 
images, since each CT scan may have different scales. In this work, 512x512 size has been used and the 
package that is used in this phase is the pillow (PIL) imaging python package version 5.4.1. Figure 7 shows 
the filtering phase results in the second box where the size and the quality of the image are changed by 
clicking on the “Filter” button. 


Load Image Filter Select & Crop Save Image Clear 


= % Jie input » project pice > n 


New folder 


v. AllFiles(*.") 


Open 


The area of The Selected Part: 


The Black area: 


Air% : 


Select & Crop Save Image 


The area of The Selected Part: 


The Black area: 


Air % : 


Figure 6. Loading a DICOM image 


2.3. Region of interest (ROD) selection 

ROI extraction can be a crucial step for many computed-aided diagnosis systems [20], [21]. The 
developed tool for this project is based on the manual segmentation; where the user of this tool can divide the 
ROI from the image by mouse clicks [22], [23]. In other word, this tool focuses on manual segmentation and 
depends on the user's knowledge and experiences in determining the ROI. To reach this phase, the user has to 
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pass through the preprocessing phase using the filter button on the developed tool as mentioned earlier. After 
showing the filtered image in the middle image box as shown in Figure 7, the user can select the ROI by 
using the “Select and crop” button. Once the user clicks this button, the CT image will be displayed in a 
separate window as Figure 8. Through this window, the user relies on mouse clicks to set the coordinates of 
the selected part as depicted in Figure 9. The selection and cropping procedures are done using the OpenCV 
python package (version 3.4.1). At the end, the tool computes the area of the ROI and the area of the air, 
which is colored black (determining the black area is done by simply inspecting the pixels’ intensities). The 
procedure used for computing the area is done using the Shapely Python package for computational geometry 
(version 1.6.4), where the area is computed using the object.area method which returns the area of the 
selected ROI. Finally, the percentage of the air is computed and displayed in a red box as Figure 9 shows. 
The tool is developed on a personal laptop with an Intel Core is CPU @ 1.80 GHz and 6 GB of RAM under 
Windows 8.1, and implemented using anaconda3-5.2.0 and python 3.6. As for the graphical user interface 
(GUD, PyQt5 was used. 


E | GUI = x 


Load Image Filter Select & Crop Save Image Clear 


The area of The Selected Part: The Black area: | Air %: 


Figure 5. The filtering phase, where the size and the quality of the loaded image are changed 


Load Image Filter Select & Crop Save Image Clear 


| The area of The Selected Part: The Black area: | Air W7 


Figure 6. The select and crop phase. A separate window showing the filtered image is used to 
select the ROI area 
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Load Image Filter Select & Crop Save Image Gear 


The area of The Selected Part: 2202.5 The Black area: 0,33098751418842226 


Figure 7. A cropped ROI with its computed air percentage 


3. RESULTS 


In this section, we present the results of our analysis of the collected dataset. The data were recorded 
in a Microsoft Excel spreadsheet and analyzed using the smart product-service system (SPSS) program 
version 16.0. Statistical significance was assessed using a two-tailed Fisher's exact test (statistical 
significance is considered for p<=0.05). The study included 50 patients (24 males and 26 females). At the 
time of the data collection, their ages ranged from 13 to 62 years with a mean age of about 29 years. Five 
patients were 16 years old or younger. 

Table 1 shows the distribution of patients over all symptoms under consideration. The right side 
nasal obstruction was the most common symptom in (58%) of the patients, followed by left side nasal 
obstruction (56%), right side headache (52%), left side headache (48%), hyposmia (42%), snoring (40%), 
right side rhinorrhea (34%) and left side rhinorrhea (30%). However, only (34%) of the patients presented 
with bilateral nasal obstruction and (36%) of patients with a headache on both sides. Table 2 shows the 
average surface area (SA) at the three different levels on both sides. At L1, which represents the SA between 
the nasal septum and the head of inferior turbinate, the average SA at the right side was 0.306 mm2 and, at 
the left side, it was 0.278 mm? Table 3 shows the patients nasal obstruction status compared to SA at L1. 


Tabel 1. The distribution of patients over all symptoms under consideration 


Symptoms Number of patients Percentage 
Bilateral Nasal obstruction 17 34% 
Right side nasal obstruction alone 12 24% 
Left side nasal obstruction alone 11 22% 
Headache on both sides 18 36% 
Right side headache 8 16% 
Left side headache 6 12% 
Rhinorrhea 12 24% 
Right side rhinorrhea 3 10% 
Left side rhinorrhea 3 6% 
hyposmia 21 42% 
Snoring 20 40% 


Tabel 2. The average SA at the three different levels on both sides 


Level Average distance 
L1 on right side 0.306 
L1 on left side 0.278 
L2 on right side 0.337 
L2 on left side 0.304 
L3 on right side 0.356 
L3 on left side 0.340 
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Tabel 3. Patients nasal obstruction status compared to SA at L1 


Avg distance at Nasal obs at Asymptomatic at # of patients Nasal obs at left | Asymptomatic # of patients 
L1 right right at left 
Below avg 16 5 21 19 8 27 
Above avg 13 16 29 9 14 23 
# of patients 29 21 50 28 22 50 


On the right side, out of 50 patients, a total of 29 patients had symptomatic nasal obstruction and 21 
were asymptomatic. Among those who were symptomatic, 16 patients had SA below the average and 13 
patients above the average while, for those who were asymptomatic, 16 patients had SA higher than the 
average and 5 patients had SA lower than the average. The difference between the two groups was 
statistically significant (p=0.042). On the left side, out of 50 patients, a total of 28 patients had symptomatic 
nasal obstruction and 22 were asymptomatic. Among those who were symptomatic, 19 patients had SA 
below the average and 9 patients above the average while, for those who were asymptomatic, 14 patients had 
SA higher than the average and 8 patients had SA lower than the average. The difference between the two 
groups was also statistically significant (p=0.0448). No other statistically significant differences were found 
when we comparing the other symptoms with the SA in L1. 

L2 represents the distance between the nasal septum and the head of middle turbinate. The average 
SA at the right side was 0.337 mm? and at the left side was 0.304 mm? Table 4 shows the patients headache 
status compared to SA at L2. Both the right and left side SA had a statistically significant correlation with the 
headache with p values 0.021 and 0.024, respectively. No other significant correlations with other symptoms 
were found at that level.Regarding L3, there was no statistically significant correlation between right or left 
SA and any of the symptoms investigated in the study. 


Tabel 4. Patients’ headache status compared to SA at L2 


Avg distance at Headache at Asymptomatic at # of patients Headache at left © Asymptomatic # of patients 
L2 right right at left 
Below avg 19 9 28 18 11 29 
Above avg 7 15 22 6 15 21 
# of patients 26 24 50 24 26 50 


4. DISCUSSION 


CT is now considered among the best diagnostic options for evaluation and assessment of paranasal 
sinuses, nasal cavity and for demonstrating various sinonasal diseases. Moreover, it is commonly used to 
monitor medical therapy for nonsurgical patients [24]. However, despite the generalized use of CT, it’s 
accuracy in diagnosing chronic rhinosinusitis is still less evident [25]. Overall, when combined with an 
endoscopic exam of the nose, it provides sufficient data that is needed for diagnosing sinuses problems 
problems [26], [27]. According to [28], for effective demonstration and visualization of nasal anatomy and 
paranasal sinuses, both endoscopy and CT should be combined. In conclusion, CT can serve as an anatomic 
map for the surgeon and it is mandatory for assessment of paranasal sinuses [26]. In our study, we present a 
CT scan as a radiological tool for evaluating patients’ symptoms with a new approach that correlates the 
symptoms with the diameter of the nasal cavity at three different levels in symptomatic patients who did not 
have changes in their CT images. Until now, there are no published articles that correlate the surface area of 
the nasal cavity using a CT scan that showed no evidence of rhinosinusitis with patients’ sinus symptoms. In 
this paper, we concluded that the patient may suffer from nasal symptoms “facial pain and nasal obstruction” 
event in case of a normal CT scan due to narrow airway passage in the nasal cavity at different levels. 

This is a descriptive clinical study carried out on 50 symptomatic sinus diseased patients who 
presented with different nasal symptoms where the nasal obstruction was the most common symptom (80% 
of the patients), followed by headache (64%), hyposmia (42%), rhinorrhea (40%) and snoring (40%). In the 
present study, the patients’ ages ranged between 13 years and 62 years. Most patients were in the third or 
fourth decades of their life with a male: female ratio of 1:1. All of them had undergone CT imaging of 
paranasal sinuses to correlate these symptoms with the different diameters of the nasal cavity at three 
different levels. Bhattacharyya et al [29] used the sinonasal outcome test 20 (SNOT-20) to compare patient- 
based symptoms with sinus CT imaging in patients referred for sinus CT. They found no correlation between 
the CT imaging and total SNOT-20 score on the facial pain or pressure question. Moreover, Stewart at al. 
[30] confirmed that there is no correlation between the radiological findings on CT scan and the symptoms 
score for rhinosinusitis. Finally, both Kenny et al. [31] and Wabnitz et al. [32] concluded that the CT scan 
findings do not account for much of the variability of the symptoms in chronic rhinosinusitis, despite the 
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statistically significant correlations that they found among symptoms and radiographic findings due to low 
coefficient of determination in their results. 

In a study done by Flinn et al. [28], the authors found that 27% of patients with no symptoms of 
rhinosinusitis had sinus opacification for those patients who referred for a head CT scan. In contrast to this, 
Calhoun et al. [33] concluded that the incidence of incidental sinus abnormalities “sinuses opacification on 
CT scan” is significantly greater in patients with a history of sinus symptoms compared to those without 
symptoms. In our study, despite the absence of opacification on sinuses, both nasal obstruction and facial 
pain were statistically significantly correlated with narrow nasal cavity passage at level one and level two 
respectively on CT scan, and this may explain the theory of absent correlation between the radiological 
finding on CT scan and the symptoms score for rhinosinusitis for [30] and [31]. We also confirmed the lack 
of correlation between snoring, rhinorrhea, and hyposmia with any of the nasal cavity surface area at the 
three different levels on CT scan. In another word, nasal obstruction is a symptom related to the narrow 
passage of air at L1 (head of inferior turbinate) while the facial pain due to narrow passage at L2 (middle 
turbinate) in CT scan patients. 


5. CONCLUSION 

In this study, we created our own dataset consists of 50 patients who visited the ENT clinics at 
KAUH. Those patients suffer from different sinonasal symptoms. Collecting the dataset is done through 
different phases and procedures based on medical image processing techniques. In addition to creating the 
dataset, we developed a tool using the python language and PyQt5 for the User Interface. This tool reads the 
DICOM images and computes the air percentage in the nasal cavity. We employ image processing techniques 
to correlate the narrowing of nasal surface area with sinonasal symptoms. In our study, we find that image 
processing techniques could be very useful tools to predict the severity of common nasal symptoms and they 
can be used to suggest treatment and to follow up on the case progression. 
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